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ABSTRACT 

Optically similar elliptical galaxies have an enormous range of X-ray luminosities. We show that 
this range can be attributed to large variations in the dark halo mass M V i r determined from X-ray 
observations. The K-band luminosity of ellipticals varies with virial mass, Lk oc 

Kir , but 

with considerable scatter, probably due to the stochastic incidence of massive satellite galaxies that 
merge by dynamical friction to form group-centered ellipticals. Both the observed X-ray luminosity 
L x oc M^ir and L x /Lk oc M^f are sufficiently sensitive to the virial mass to explain the wide variation 
observed in L x among galaxies of similar Lk- The central galaxy supernova energy per particle of 
diffuse gas increases dramatically with decreasing virial mass and elliptical galaxies with the lowest 
X-ray luminosities (and M V i r ) are easily explained by supernova-driven outflows. 

Subject headings: X-rays: galaxies - galaxies: clusters: general - X-rays: galaxies: clusters - galaxies: 
cooling flows 



1. INTRODUCTION AND BACKGROUND 

In spite of their very regular optical properties, ellip- 
tical galaxies with similar optical luminosities have X- 
ray luminosities that range over more than two orders 
of magnitude. Many attempts to explain this huge scat- 
ter or to find observable properties that correlate with 
residuals in this scatter have been inconclusive (e.g. Es- 
kridge, Fabbiano & Kim 1995a, b,c; Ellis & O 'Sullivan 
2006; Mathews & Brighenti 2003 for a review). 

Figure I shows a plot of X-ray and K-band luminosi- 
ties for elliptical galaxies taken from the catalog of Ellis 
& O 'Sullivan (2006). Most of the X-ray emission is from 
diffuse hot gas although stellar emission from low mass 
X-ray binaries (LMXB) contributes in galaxies of the 
lowest L x . While the lower envelope of the distribution 
is well-defined by the locus for stellar LMXB emission, 
X-ray images of these \ow-L x elliptical galaxies also show 
definite evidence of diffuse gas emission (e.g. Sarazin et 
al. 2001). There is no apparent concentration of gas- free 
ellipticals along the LMXB locus, confirming the general 
prevalence of diffuse emission. Figure 1 includes only a 
few ram-stripped, low-L x elliptical galaxies orbiting in 
massive clusters that have nevertheless retained a small 
mass of - IkeV gas (e.g. Sun et al. 2004; 2005). Our 
interest here is not these strongly ram-stripped galaxies, 
but similar galaxies that are isolated or central galaxies 
in groups. 

Using an early data set based on Einstein, ROSAT 
and ASCA observations, Mathews & Brighenti (1998a) 
showed that the radius of half- X-ray emission r e ff x de- 
creases systematically with decreasing L x in Figure 1. It 
was demonstrated later that this correlation could be un- 
derstood as a variation of the dark halo mass (Mathews 
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& Brighenti 1998b). In this paper we extend this idea 
using more recent data to show that the scatter in L x is 
intimately related to stochastic and systematic variations 
of the virial mass M OT > of the dark halos that surround 
them. Galaxies with smaller M OTr necessarily have fewer 
baryons and lower L x and we use simple galaxy plus halo 
models to show that this effect alone can explain at most 
about half of the total observed range in \ogL x . How- 
ever, galaxies/groups with the lowest L x require AGN or 
supernova-driven outflows. 

2. CORRELATIONS AMONG L x , L K AND M VIR 

The virial mass of dark host halos around elliptical 
galaxies can be directly determined from radial temper- 
ature and density profiles in the hot gas, assuming hy- 
drostatic equilibrium. The total mass based on Chandra 
and XMM images can be attributed to a combination 
of stellar and dark mass profiles. Chandra observations 
are particularly valuable since LMXB sources can be re- 
solved and removed. The galaxy sample we consider here 
with known viral masses is shown in Figure I with filled 
symbols: squares are from Humphrey et al (2006), cir- 
cles are from Gastaldello et al. (2006) and the hexagons 
are X-ray luminous galaxies tabulated in Mathews et al. 
(2005). 

Figure 2a shows that M OT > for the sample galaxies in- 
creases with L x according to the correlation 

\ogM mr = -(4.20 ±0.92) + (0.42 ±0.02) log L x (I) 

(with correlation coefficient 0.75) indicated with a dashed 
line, i.e. L x oc M^. This and subsequent correlations 
are formal least square fits, assuming a mean typical er- 
ror A log = 0.090 which is small compared to the 
intrinsic spread in the data. The quality of our least 
squares fit in Figure 2a is limited by the cosmic spread 
in M V i r , but it is consistent with a slightly steeper varia- 
tion M v i r oc L x in massive clusters found by Reiprich 
& Bohringer (2002). 

Figure 2b indicates a very wide variation in the virial 
mass with K-band luminosity of the central galaxy Lk, 
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superimposed on a positive M vir -LK correlation, 
log AW = -(1.75 ±1.04) + (1.34 ± 0.09) log Lk, (2) 

(dashed line in Fig 2b) with correlation coefficient 0.58. 

The mean optical luminosity of group and cluster- 
centered elliptical galaxies is known to increase, then 
saturate as the virial mass of the cluster increases. Halo 
masses around elliptical galaxies in small groups can be 
determined from the SDSS galaxy-mass correlation func- 
tion (Cooray and Milosavljevic 2005), from weak galaxy- 
galaxy lensing (Yang et al. 2003), and from the X-ray 
temperature of the gas (Lin & Mohr 2004). The Lk of 
the central galaxy varies as M ^ 5 for M vir < 10 12 M Q 
and as M^f for M vir > 10 14S M Q . The intermediate 
variation Lk <x M°-J 5 in equation (2) is consistent with 
these results for the range in M vir that we consider. 

The scatter in Figure 2b arises because of an intrin- 
sic (log-normal) cosmic variance of M V i r among galaxies 
of similar Lk- The central galaxy forms from a rather 
small number of luminous satellite galaxies that merge at 
the halo center because of dynamical friction (e.g. Fal- 
tenbacher & Mathews 2005). Since dynamical friction 
varies as the square of the ratio of the galaxy subhalo 
mass to the mass of the host halo, only the most luminous 
group satellite galaxies merge at the center. Therefore, 
because of the relative rarity and stochastic incidence of 
luminous satellite galaxies, halos of the same mass can 
have a range of Lk for the central elliptical (e.g. Cooray 
2006). Conversely, halos of different M„j r can contain 
central galaxies with similar Lk, as seen in Figure 2b. 

Figure 2c shows that the mean gas temperature (kT) 
(estimated from T(r) profiles for the Humphrey et al 
2006 data) is tightly correlated with the virial mass. The 
dashed line describing this correlation is given by 

log M OTr = (13.5 ± 0.02) + (2.07 ± 0.09) log(fcT). (3) 

This kT-M V i r correlation is similar to those found in 
clusters with more massive halos (Arnaud et al. 2005; 
Vikhlinin et al. 2006). 

3. RELATION OF L x TO BARYON CLOSURE 

The physical origin of the vertical deviation of L x in 
Figure 1 is best characterized by the logarithmic distance 
D(\ogL x ) of each galaxy from the (dash-dotted) locus 
of baryonic closure, the approximate L x expected if the 
most massive dark halos for each Lk are filled with gas 
(and stars) containing the maximum baryonic mass al- 
lowed by the cosmic baryon ratio /(, = f2f,/f2 TO = 0.176 
observed by the Wilkinson Microwave Anisotropy Probe 
(WMAP; Spergel et al. 2006). If a galaxy/group is bary- 
onically closed, no gas has been ejected from the confin- 
ing halo by AGN or supernova energy. Since the total 
mass of diffuse gas can be reasonably well estimated in 
several of the most luminous groups in the sample, we 
know that at least some galaxies/groups with the largest 
L x are quite close to the closure locus (Mathews et al. 
2005), 

(\ogL x ) c i osure w 18.361 + 2.18 log L K , 

shown by the dash-dotted line in Figure 1. The deviation 
in L x of a given galaxy below the closure locus is 

D(\0gL x ) = (log L x ) closure ~ \ogL x . (4) 



Figure 2d shows an anticorrelation between D(\ogL x ) 
and M V i r 

log M mr = (14.01 ± 0.03) - (0.34 ± 0.03)-D(logI, x ), (5) 

(dashed line in Fig. 2d) , confirming that the virial mass 
is key to understanding the wide scatter in L x among 
elliptical galaxies in Figure 1 . Eliminating M V i r between 
equations (1) and (2), we derive the approximate direc- 
tion in the L x -Lk plane along which M V i r decreases most 
rapidly, shown by the short-dashed arrow in Fig. 1. 

4. ISOTHERMAL MODELS 

The closure locus in Figure 1 - the approximate max- 
imum M v i r and L x for galaxy/groups for each Lk - 
can be related to simple models for the hot gas. We 
assume that the gas is isothermal with {kT) given by 
equation (3) and in hydrostatic equilibrium in the po- 
tentials of the dark halo and central galaxy. The NFW 
halo is characterized by its virial mass M V i r (using a 
mean concentration c{M V i r ) from Bullock et al. 2001). 
The virial radius is defined by an overdensity A = 104 
above the critical density p c — 9.24 x 10~ 30 gm cm~ 3 . 
The central galaxy is assumed to have a de Vaucouleurs 
mass distribution described by Lk and effective radius 
r eff = 4.16(M»/10 n M Q ) - 56 kpc (Boylan-Kolchin et al. 
2005). The B-band stellar mass to light ratio is T B = 
M*/L B = 7.5 and we adopt L k /L K q = 4:.15L B /L BQ 
appropriate for an old stellar population in elliptical 
galaxies (Ellis & O'Sullivan 2006). Solutions of these 
isothermal, hydrostatic models with fb = 0.176 (and in- 
cluding gas self-gravity) give the maximum L x within 
r V i r for any given M„j r and Lk- For each Lk we increase 
M vir until L x of the model exceeds all observed values 
and eventually recover the approximate closure locus de- 
scribed by Mathews et al. (2005) (dash-dotted line in 
Figure 1), given by L x = 2.19 x lO^L^/lO 11 L K q) 2A8 
erg s _1 . The closure locus empirically defines an upper 
limit to the L x -Lk distribution and is consistent with the 
large baryon fractions « f b determined for well-observed 
galaxy/groups having the largest L x /Lk- Since L x of 
the models extend to r vir while the observed L x do not, 
even galaxies with the largest L x in Figure 1 are expected 
to lie slightly below the closure locus. 

In principle, galaxies that lie significantly below the 
closure line in Figure 1 could also be baryonically closed, 
i.e. with fixed ft = 0.176, L x is expected to decrease with 
M v i r . Although we view this possibility as highly un- 
likely, it is interesting to explore how much of the spread 
in Figure 1 can be understood with this extreme assump- 
tion of universal baryon closure. Consider simple models 
in which M v i r decreases below the closure locus in Figure 
1, holding Lk of the central galaxy fixed and ignoring the 
stellar mass of satellite galaxies. For central galaxies of 
mass M* = T B L B = (7.5/4.15)L^ M the fraction of 
baryons remaining in diffuse gas is 

fa = M gas /M vir = f b - 1.8L K /M mr . (6) 
L x declines with M v i r faster than M gas since stars (with 
constant Lk) consume an increasingly larger fraction of 
the baryons in smaller halos. In this simple model, where 
mass loss from evolving stars is neglected, the mass of 
diffuse gas (and L x ) completely vanishes (f g — 0) for 
halo masses less than M V i r ^ — IO.SLk ■ 

Isothermal models with logi^- = 11.5 and 
log(M mr /M a ) = 13.5, 13 and 12.6 (models a, b, and 
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c respectively) are listed in Table 1 and plotted in Fig- 
ure 1. From equation (6) model c is close to the f g = 
limit. The X-ray luminosities of these models extend over 
nearly half the total spread in Figure 1 . The radius r e f f x 
in Table 1 that includes half of the total projected L x de- 
creases systematically with L x , qualitatively confirming 
the correlation between the observed X-ray sizes and L x 
found by Mathews & Brighenti (1998a). Also listed in 
Table 1 are the thermal and binding energies per gas par- 
ticle, eth and eund, as well as the corresponding energy 
released by supernovae of Type II, e/j, calculated assum- 
ing a Salpeter IMF, and Type la, e/ a , estimated from the 
total iron mass produced by an old stellar population, 
Mpe/M* w 0.0006 (De Grandi et al. 2004), assuming 
0.7 Mq in iron per supernova. 

As M vir decreases, the supernova energy in the cen- 
tral galaxy varies from relative insignificance for model 
a to dominance in model c where en + e/ a 3> etind, the 
gravitational binding energy per gas particle. The global 
energetics in Table 1 strongly suggest that the remain- 
ing scatter in L x in Figure 1 below model c can easily 
be explained by supernova-driven gaseous outflows from 
low M v i r halos. It is more likely that outflows reduce 
L x for all galaxies in Figure 1, and that outflows are dis- 
proportionally effective as M V i r decreases with increasing 
D(\ogL x ) (Fig. 2d). 

Another feature in the galaxy distribution in Figure 1 
that remains to be understood is the absence of galaxies 
near the closure locus but with logL^ < 11.2. If the 
halo mass decreases in the direction of the short-dashed 
arrow, it is plausible that the absence of galaxies in this 
region can also be understood in terms of supernova- 
driven outflows combined with a modest gas supply from 
evolving stars, maintaining L x close to the LMXB lower 
limit. 

Finally, we note another observable feature that corre- 
lates with M v i r : the shape of the radial gas temperature 
profile T{r). Massive galaxy groups with large L x have 
a temperature maximum displaced from the center, but 
in low-Lz ellipticals the temperature peaks at the cen- 
ter (e.g. Humphrey et al. 2004; Humphrey et al. 2006; 
O' Sullivan & Ponman 2004). The origin of this difference 
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in thermal profiles can be explained by the ratio of the 
mass of the galaxy to its host halo (Mathews & Brighenti 
2003). When a group-centered elliptical resides in a more 
massive halo, the virial (and gas) temperature through- 
out most of the halo exceeds that of the central galaxy, 
resulting in a positive temperature gradient in the inner 
regions. However, when a similar elliptical has a halo 
of lower mass, the gas temperature peaks at the center 
where it approaches the virial temperature of the stel- 
lar bulge. This observed pattern is consistent with the 
expected variation of L x with M v % r in Figures 1 and 2a. 

5. CONCLUSIONS 
Our main conclusions are: 

(1) The large variation of X-ray luminosity among group- 
centered elliptical galaxies of similar optical luminosity 
can be ascribed to a large cosmic variation in the sur- 
rounding dark halo mass. 

(2) For galaxies with similar Lk a significant fraction 
of the observed correlation between L x and M v ^ r can 
be expected simply because smaller halos contain fewer 
baryons, less hot gas, and lower L x . However, with de- 
creasing M V i r /LK supernova (or AGN) driven outflows 
become more likely, reducing L x further. 

(3) M vir decreases with both L x and Lk in the L x -Lk 
plot. 

(4) The upper envelope of elliptical galaxies in the L x - 
Lk plot lies close to the locus of WMAP baryon closure 
for galaxies having the most massive dark halos at each 
L K . 

(5) Previously found correlations between the X-ray lu- 
minosity of elliptical galaxies (having similar Lk) and 
their X-ray sizes (Mathews & Brighenti 1998a) or tem- 
perature profiles (Humphrey et al. 2006) are also consis- 
tent with a systematic decrease in halo mass for smaller 
L x . 
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TABLE 1 

ISOTHERMAL MODELS WITH L K = 10 



model 


a 


b 


c 


log M vir (M Q ) 


13.5 


13.0 


12.6 


log kT (keV) 


1.0 


0.6 


0.4 


log M gaa (M @ ) 


12.66 


12.07 


11.11 


log (erg/s) 


43.43 


42.77 


42.09 


r e ffx ( k Pc) 


151 


104 


0.9 


e th (kcV/part) 


1.5 


0.9 


0.6 


ebind (kcV/part) 


2.3 


1.1 


0.6 


en (keV/part) 


0.4 


1.5 


14. 


e /a (keV/part) 


0.03 
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Fig. 1. — K-band and X-ray luminosities of X-ray-detccted elliptical galaxies (T < —4) from the catalog of Ellis & O'Sullivan (2006). 
Filled squares, circles and hexagons are our sample galaxies for which the virial mass has been determined except RXJ146 (the open 
hexagon) which is useful to establish the upper envelope of the data. The lower long-dashed line shows the approximate emission from 
low mass X-ray binary stars (Kim & Fabbiano 2004). The upper dash-dotted line is the approximate locus of baryon closure (Mathews et 
al. 2005). The short-dashed arrow shows the direction in which M v i r decreases most rapidly. Labels a, b and c show locations of several 
baryonically closed, isothermal models listed in Table 1. 



6 



SO 
O 



44 
43 
42 
41 



40 
11.8 
"J 11.6 



Clfl 
O 



A 
H 

V 
O 



11.4 
11.2 
11 

0.5 b- 



-0.5 



X 

J 1 




: □ 


o 


L □ 

: n 


Q 





i i i — r~| — n — r~r~| — r~i — r~i — prn — pn — p 



a 



_ □ CP. '° 

KM 



I I I I I I I I I I > I 



o.9g o 



□ □ 







□ 



o 







I I I l x l I I I I I I I I I I I I 



□ 



'□ 



t ...... ^.....y 



d 



1% ° 



J I I I ' I I L_ 



I I I I I I 



12.5 13 13.5 14 14.5 

log M vir (M s ) 



vir v " 0' 

— 1 ) of the sample galaxies with the virial mass (in Mq) of their host halos. 



Fig. 2. — (a) Variation of X-ray luminosity L x (in erg s~ 
(b) Variation of the central galaxy K-band luminosity (in Lkq) with virial mass, (c) Variation of mean gas temperature (in keV) with 
virial mass, (d) Variation of the logarithmic deviation from the locus of baryon closure in Figure 1 with the halo virial mass. Note that 
D(logL x ) increases downward, as in Figure 1. Dashed lines show correlations discussed in the text. Note that all correlations are apparent 
in published data from Humphrey ct al (2006) (open squares) and Mathews ct al. (2005) (open hexagons), while the data of Gastaldello ct 
al. (2006) (open circles) fall in between. 



